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ABSTRACT 

In a search for the signature of turbulence in the diffuse interstellar medium in gas density 
distributions, we determined the probability distribution functions (PDFs) of the average vol- 
ume densities of the diffuse gas. The densities were derived from dispersion measures and 
Hi column densities towards pulsars and stars at known distances. The PDFs of the average 
densities of the diffuse ionized gas (DIG) and the diffuse atomic gas are close to lognormal, 
especially when lines of sight at |6| < 5° and |6| ^ 5° are considered separately. The PDF 
of («.Hi ) at high |6| is twice as wide as that at low The width of the PDF of the DIG is 
about 30 per cent smaller than that of the warm H i at the same latitudes. The results reported 
here provide strong support for the existence of a lognormal density PDF in the diffuse ISM, 
consistent with a turbulent origin of density structure in the diffuse gas. 
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1 INTRODUCTION 

Simulations of the interstellar medium (ISM) have shown that, if 
isothermal turbulence is shaping the structure of the medium, the 
density distribution becomes lognormal fElmegreen & Scalo|2004| 
and references therein). However, Madsen et al. ( 2006J found large 
temperature differences between lines of sight through the warm 
ionized medium, which seem inconsistent with an isothermal gas. 
Also in the MHD simulations of |de Avillez & BreitschwerdtH2005^ 
and |Wada & Norman|p007"l l the ISM became non-isothermal. The 
latter authors showed that for a large enough volume, and for a long 
enough simulation run, the physical processes causing the density 
variations in the ISM in a galactic disc can be regarded as ran- 
dom and independent events. Therefore, the PDF of log(density) 
becomes Gaussian and the density PDF lognormal, although the 
medium is not isothermal. The medium is inhomogeneous on a lo- 
cal scale, but in a quasi-steady state on a global scale. 

The shape of the gas density PDF can be an important com- 
ponent in theories of star formation. |ElmegreenH2002) showed that 
a Schmidt-type power-law relation between the star formation rate 
per unit area and the gas surface density can be deduced if the den- 
sity PDF is lognormal and star formation occurs above a threshold 
density: the dispersion of the lognormal density PDF is a key pa- 
rameter in the model. However, from the existing simulations it is 
not yet clear whether a lognormal density PDF in galaxies is uni- 
versal and which factors determine their shape jWada & Norman] 
[2007l >. 

Little observational evidence exists to test the results of the 



simulations. |Wada et al.| j2000| l showed that the luminosity func- 
tion of the HI column density in the Large Magellanic Cloud is 
lognormal. Recently, pill et al.|p007||2008 !) derived a lognormal 
distribution of the emission measures perpendicular to the Galac- 
tic plane of the DIG in the Milky Way, observed in the Wisconsin 
Hq Mapper survey jHaffner et aL|2003) above Galactic latitudes 
of 10°. They note that emission measures towards classical HII re- 
gions also fit a lognormal distribution, with different parameters. 
Furthermore, Tabatabaeij l|2008| found a lognormal distribution of 
the emission measures derived from an extinction-corrected H-a 
map of the nearby galaxy M33 ( jTabatabaei et al.||2007^ . [Gaustad| 
|& Van Buren| ( [T993] l plotted the distribution of the local volume 
density of dust near stars within 400 pc of the Sun, which is also 
consistent with a lognormal distribution. In this Letter we present 
the PDFs of average volume densities of the DIG and the diffuse 
atomic gas in the Milky Way, and show that they are consistent 
with lognormal distributions as well. 



2 BASICS AND DATA 

We investigated the density PDFs of the DIG and of diffuse atomic 
hydrogen gas (H i) in the solar neighbourhood. 

Various volume densities of the DIG can be obtained from the 
dispersion measure DM (in cm"'' pc) and emission measure EM 
(in cm"*^ pc) towards a pulsar at known distance D (in pc) from 
the relations: 
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where Uc (l) is the local electron density at point / along the line-of- 
sight (LOS), {ric) and (n^) are averages along D and Fd is the 
fraction of the line of sight in clouds of average density A^c (see 
fig. 1 in|Berkhuijsen et al. 20061. The final equality in Eq.|2]is only 
valid when the average density of every cloud along the LOS is 
the same: then ( ) = (n^ ) Fd = F^. Thus A^c and Fd are 
crude approximations of the true average cloud density and filling 
factor along a LOS, but given the large number of different LOS in 
our sample their mean and dispersion are reasonable estimators of 
iVc andFd. 

Combining Equations and ^ we hav^ 
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The line-of-sight filling factor Fd approximates the volume filling 
factor Fv if there are several clouds along the line of sight. As this 
will generally be the case, we take Fv = Fd. 

We used the densities derived from two pulsar samples that 
were originally selected for studies of the volume filling factor of 
the DIG. 

(i) 34 pulsars at observed distances known to better than 50 per 
cent, collected by Berkhuijsen & Muller| ^2008) . The pulsar dis- 
tances are in the range 0.1 < D < 9.5 kpc, with a mean distance 
of 2.4 kpc and a standard deviation of 2.9 kpc (the spread is large 
because 21 pulsars lie at D < 2 kpc). 

(ii) 157 pulsars with distances obtained from DM and the model 
of the distribution of free electrons in the MW of Cordes & Laziol 
( |2002^ , collected by [Berkhuijsen et al] j2006[ l. These pulsars lie 
in the range 0.1 < D < 6 kpc, with mean distance L7kpc and 
standard deviation 1.0 kpc. 

Apart from 6 pulsars in the small sample, all pulsars are located at 
1 6| 5° in order to ensure that the lines of sight towards the pulsars 
probe the DIG and not denser H ii regions. The dispersion measures 
were taken from the catalogue of [Manchester et aT] ^2005l l. The 
emission measures in the direction of the pulsars were obtained 
from Hq surveys ( [Haffner et al. 2003 ; Finkbeii iere't al.|2002) and 
corrected for extinction piplas^fe Savagejl994 b, Dickin son et al.| 
|2003j > as well as for Ha emission originating beyond the pulsars. 
We refer to the work of [ Berkhuijse n et al.|p006| hereafter called 
BMM) and 'Berkhuijsen & Muller| ^2008| for further details[^ 

Diplas & Savage ( 1994b) studied the scale height of the dif- 
fuse dust and the diffuse atomic gas using 393 stars in the Galaxy. 
We calculated the average Hi volume density, { tihi ) (in cm~'^), 
from the column density A'^(Hi) (in cm~^), corrected for contri- 
butions from the star, and the distance to the star as given in Table 
I of |Diplas & Savage| l fT994a) : 
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where nni{l) (in cm '') is the local volume density at distance / 



Figure 1. Probability distribution functions of mean density in clouds, A^c 
(a and c) and volume filling factors, Fv (b and d). The dashed lines are 
the lognormal fits to the histograms; the fit parameters are given in Table[T] 
Top: small sample (34 pulsars); bottom: BMM sample (157 pulsars). 



along the line of sight. We removed 18 stars from the sample with 
denser clouds in their lines of sight indicated in Fig. 9 of |Diplas"&| 
|Savage|(l994b^ , leaving the data towards 375 stars for analysis. The 
stars in this sample have distances in the range 0.1 < D < 11 kpc, 
with a mean distance of 2.2 kpc and standard deviation 1.7 kpc.|^ 



3 RESULTS 

3.1 Density PDFs of the DIG 

In Fig. [T^ we present the probability distribution function (PDF) 
of the mean density in clouds, TVc, for the sample of 34 pul- 
sars. As the volume filling factor is (anti-) correlated with A^'c (see 
BMM), we show the PDF of Fv for the same sample in Fig. lb. In 
log space, both PDFs are consistent with a Gaussian distribution, 
which is equivalent to a lognormal distribution in linear space. The 
PDFs have about the same dispersion, a (see Table [TJ. The posi- 
tions of the maxima, /i = log(density of maximum), correspond to 
A^c = 0.19 ± 0.02 cm~^ and Fv = 0.078 ± 0.006, which repre- 
sent the centre of gravity in the Fv - A'^c plot in fig. 6 of Berkhuijsen] 
|& Mtil ler (2008). This sample is rather small, with low counts A'^; 
in the histogram bins and (probably) overestimated Poisson errors 
5i — \fNi leading to rather small reduced-x^ statistics (Table[T]l. 
Therefore we also calculated the PDFs of Nc and Fv for the much 
larger sample of BMM, which are shown in Figs. lc,d. Both are 
well fitted by Gaussians of widths that are nearly identical to those 
of the small sample (see Table 1). The positions of the maxima are 
at A^c = 0.19 ± 0.01 cm-^ and Fv = 0.115 ± 0.005, corre- 
sponding to the centre of gravity in the Fv- A'^c plot of BMM (their 
fig. 11). The good agreement between the PDFs of the two samples 
indicates that the statistical results on Fv and A'c of BMM are not 
influenced by the model distances and statistical absorption correc- 
tions that they used. 

In Fig. [2] we present the PDFs of the average densities { Tie ) 
and { Tie ) for both samples, all of which are well described by a 
lognormal distribution. The dispersion in ( Tie ) is smaller than the 



^ Note that we write Fd , Fv and A^c where [Berkhuijsen et al.|j2"006} used 
/d, /v andrfc. 

^ The data used in our analysis are available, on request, from EMB. 



^ We do not consider the PDFs of the column densities DM, EM and 
A'^(H l), because they are influenced by the distributions of the distances to 
the pulsars and stars in the samples. 
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Table 1. Lognormal fits to the PDFs of the DIG (Figs.fTlandB. The fitted function isY = (v^cr)-l exp[- (logio ^ - m)^/2o-2]. 







Position of maximum 


Dispersion 




Sample 


X 




X 


fj 




This work 


Nc (cm-3) 


-0.71 ± 0.03 


0.19 ±0.02 


0.32 ±0.03 


0.4 


N = 34 




-1.11 ± 0.03 


0.078 ± 0.006 


0.31 ± 0.03 


0.8 




(ric ) (cm~^) 


-1.83 ±0.02 


0.015 ±0.001 


0.22 ±0.02 


0.3 




(n2) (cm-6) 


-2.48 ±0.02 


0.0033 ± 0.0002 


0.29 ±0.02 


0.4 


BMM 


Afc(cm-3) 


-0.72 ±0.02 


0.19 ±0.01 


0.29 ±0.02 


1.7 


N = 157 




-0.94 ±0.02 


0.115 ±0.005 


0.33 ±0.02 


0.9 




(ric > (cm~^) 


-1.68 ±0.01 


0.021 ± 0.001 


0.12 ±0.01 


1.2 




(n2) (cm-6) 


-2.41 ± 0.03 


0.0039 ± 0.0003 


0.34 ± 0.03 


0.8 



Y is the fraction of sightlines in each bin divided by the logarithmic bin-width d(log]^g X). is the reduced chi-squared goodness of fit parameter, with the 
error in each bin Si estimated as Si = ^/Nl and for (number of bins— 2) degrees of freedom. 
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Figure 2. Probability distribution functions of ( ric ) = DM/D(aandc) 
and {riD = EM/D (b and d) for the small sample (top) and the BMM 
sample (bottom). The dashed lines are the lognormal fits to the histograms; 
the fit parameters are given in Table[T] 



0.4 
0.3 

J 

SO.2 

T 

'o.i 



0.4 
0.3 

>« 
J 

io.2 

T 

"0.1 



-1 
log( <n[HI]> ) 



-1 
log( <n[HI]> ) 





a 


b<5" 




b 










A 




2 -1 





-1 




1 


warm 


c 


b^5• 




d 




v ■ 









0.4 
0.3 

0.2 ■ 

I 

0.1 " 


0.4 
0.3 
0.2 \ 
0.1 " 



Figure 3. Probability distribution functions of ( nn i ) = A''(H l)/D: a) 
Full sample of 375 stars jDiplas & Savage|1994a) , of which are b) 213 
stars at \h\ < 5° and d) 162 stars at |6| > 5°; c) 140 stars seen through 
warm H I. The dashed lines are the lognormal fits to the histograms; the fit 
parameters are given in Table|2] 



dispersions in i<V and A^c due to their (anti-) correlation: and 
A'c are not independent random variables. Note that the dispersion 
of the PDF of ( Tie ) of the BMM sample is about half that of the 
small sample. As BMM used distances to the pulsars derived from 
the NE2001 model of .Cordes & Lazio| ^002; !, (n^) = DM/D 
returns the densities of the model. The small dispersion reflects the 
fact that the model is much smoother than the density variations 
in the real ISM measured for the small sample. The dispersion in 
( rie ) is larger than that of {ric) as the intrinsic spread in EM is 
much larger than in DM (see BMM; |Berkhuijsen & Miiller|2008^ 
while the distances used to calculate { ) and ( n,e ) are the same. 

It is interesting to compare our data with the results of the 
magneto-hydrodynamic simulations of the ISM in the solar neigh- 
bourhood made by |de Avillez & Breitschwerdt|pO()5} . Their fig. 
7 shows the density PDFs of five temperature regimes that devel- 
oped after about 400 Myr. The curve for 8000 < Te < 16000 K, 
which is most applicable to the DIG, closely resembles a lognor- 
mal with maximum at log(nH) = —0.75 and dispersion 0.52. The 
lognormal distribution extends over a much larger density range 
(at least —2.5 < log(nH) < 1-2) than our observations of A^'c 
(—2 < log(7Vc) < in Fig. 1). The density of the maximum of 
0.18 cm"^ agrees well with that of iVc = 0.19 ± 0.02 cm"^ de- 



rived by us (see Table 1), but the dispersion is about 70 per cent 
larger. This could be due to the larger temperature range of this 
component in'de Avillez & Breitschwerdt ( 2005 | compar ed to the 
6000-10000 K observed for the DIG i jMadsen et al.|2006| >. 

We conclude that the PDFs of the electron densities and filling 
factors in the DIG in the solar neighbourhood are lognormal as is 
expected for a turbulent ISM from numerical simulations. 



3.2 Density PDFs of diffuse H i 

In Fig. |3ji we present the PDF of the average volume density of 
Hi, (riHi ), for the full sample of 375 stars of Diplas & Savage] 
l |1994a[ >. Above log(nHi ) = — 1 the distribution is approximately 
lognormal, but there is a clear excess at lower densities reflected in 
the large reduced-x^ statistic (see Table [2|. Because low densities 
can be expected away from the Galactic plane, we calculated the 
PDFs for the latitude ranges |6| < 5° and \b\ ^ 5° separately, 
as shown in Figs.jsj) and[3|i. Both distributions have a lognormal 
shape but they are shifted with respect to each other: the maximum 
of the low-|6| sample is at { nni ) = 0.30 ± 0.02 cm""' and that of 
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Table 2. Lognormal fits to the PDFs of ( riH I > (Fig.ji}. The fitted function is Y = (v^o-)-l exp[-(logio( "H I > - m)^/2o-^]. 







Position of maximum 


Dispersion 




Area 


N 


M 


("HI > 


fj 




All stars 


375 


-0.57 ± 0.03 


0.27 ± 0.02 


0.34 ± 0.03 


6.0 


|6| < 5° 


213 


-0.52 ±0.02 


0.30 ± 0.02 


0.24 ± 0.02 


2.7 


\h\ ^ 5° 


162 


-0.80 ± 0.02 


0.16 ±0.01 


0.50 ± 0.02 


0.4 


Warm Hi 


140 


-0.98 ±0.03 


0.10 ±0.01 


0.28 ±0.03 


2.5 


\h\ < 5° 


42 


-0.89 ±0.01 


0.13 ±0.01 


0.18 ±0.01 


0.1 


\b\ ^ 5° 


98 


-1.07 ± 0.03 


0.09 ±0.01 


0.33 ±0.03 


1.0 



Y is the fraction of stars in each bin divided by the logarithmic bin-width d(logxo( I })• is reduced chi-squared goodness of fit parameter, with the 
error in each bin hi estimated as hi = vTVj and for (number of bins— 2) degrees of freedom. 



the high-j&l sample at ( nni ) = 0.16 ± 0.01 cm"^ (see Table 2). 
The latter sample clearly causes the low-density excess in Fig.[3ji. 

The dispersion of the PDF of the high-jfej sample is twice that 
of the low-jbj sample. It is not clear whether this is a real differ- 
ence or due to selection effects in the low-|6| sample. Stars at low 
Galactic latitudes can only be seen through holes between the many 
dust clouds and the low latitude sample may be biased towards low 
densities if the higher density diffuse gas ((nni) > lcm~^) is 
associated with these clouds. 

[Diplas & Savagel (199451 fig. 9) identified about 140 lines 
of sight probing the warm diffuse Hi. This sample is especially 
interesting for comparison with the DIG because the average gas 
temperatures of both components are ~ 8000K. The PDF of the 
warm H i (see Fig. |3j;) is also lognormal and peaks at { wh i ) = 
0.10 ± 0.01 cm"^ (see Table 2). All densities are < 0.3 cm"^. 
Comparison with the full sample in Fig. [3^ shows that all lines of 
sight with log ({ nni )) < —0.8 (or (nni) < 0.16 cm^^) probe 
the warm Hi. Table |2] shows that the dispersions of the PDFs of 
the warm H i are slightly smaller than for the full sample. Clearly 
the combination of warm and cool (denser) gas in the full sample 
increases the dispersion because the density range becomes larger. 



4 DISCUSSION AND CONCLUSIONS 

The results in Sect. |3] show that the average volume densities of 
the DIG and the diffuse H i within a few kpc of the Sun follow a 
lognormal distribution, as is expected if the density is the result of 
a random, nonlinear process such as turbulence. 

The dispersions of the observed PDFs vary between about 0.2 
for the DIG and for HI at |6| < 5°, and 0.5 for HI at |6i ^ 5°. Can 
we understand such differences in the frame of the simulations? 

The most remarkable difference is that the dispersions of the 
density PDFs of HI at |6| 5° are about twice those at |b| < 5° 
(see Table |2]l. If this is a real effect (see Sect. |3.2^ , a possible ex- 
planation is that the low latitude LOS cross more turbulent "cells", 
where the size of a cell is related to the decorrelation scale of the 
turbulence, than at high latitudes. Under the central limit theorem, 
one expects the density PDF to become narrower as the number 
of cells along the line-of-sight in a sample increases ( Vazquez- 
[Semadeni & Garcia|2001) l. Since the average distance to the stars 
is the same in both samples, this would imply that the average size 
of a cell is smaller at low latitudes. This would be consistent with 
the inverse dependence of the volume filling factor on mean 
density in cells/clouds found for the DIG ( [Berkhuijsen et al.|2006| 



Berkhuijse n & Muller|2008^ and diffuse dust (Gaustad & Van Bu-| 
|ren|1993| ). Higher density means smaller hence smaller clouds 
(as fewer clouds at low jfcj is unlikely). |Vazquez-Semadeni & Gar^ 
|cia| ( [200 T| used models of isothermal turbulence to investigate the 
relation between the shape of average density PDFs and the num- 
ber of turbulent cells along the line-of-sight. Their results suggest 
that the dispersion in density will narrow by a factor of 2 if the 
number of cells increases by a factor of ~ 5. 

Another interesting difference exists between the dispersions 
of the sample of warm HI at |6| ^ 5° and that of { nc ) of the DIG 
(small sample), which is at the same latitudes. The temperatures of 
the two components are similar and if the ionized and atomic gas 
are well mixed, one would expect their dispersions to be the same. 
However, the dispersion of the DIG sample, 0.22 ± 0.01, is about 
30 per cent smaller than that of the warm H i sample, 0.33 ± 0.03 
(see Tables [T] and |2](. A plausible explanation for the difference, 
which is also consistent with the higher density of the maximum in 
the diffuse { nn i ) PDF, is that low density regions are more read- 
ily ionized than higher density gas and that the average degree of 
ionization of the diffuse gas is substantially lower than 50 per cent. 
We estimate the degree of ionization to be about 14 per cent, us- 
ing the densities for { Uc ) (small sample) and ( nn i ) in Tables [T| 
and|2] consistent with the results of |Berkhuijsen et ar] l |2006[ their 
fig. 13) for a mean height above the mid-plane of about 500 pc. Al- 
ternatively, the DIG could have a higher mean temperature than the 
warm H i but with a smaller temperature range; in the simulations 
of de Avillez & Breitschwerdt (2005) high temperature gas indeed 
has a lower median density and smaller dispersion. 

Several groups have noted a link between the rms-Mach num- 
ber M and the dispersion of the gas density PDF in isothermal 
numerical simulations (e.g.lPa doan et al.|I997[|Passot & Vazquez-| 
Semadeni 1998; Ostriker et al.'|2001|l. Although the DIG is not 



isothermal ^Madsen et al.,2006j l, to a first approximation it can be 
considered so because the sound speed scales as Cs ~ y^/^ ^nd 
the observed temperature range of 6000 K < T < 10000 K cor- 
responds to only a 30% difference in Cs- Th en using the formul a 
o-fn = ln(l + /3^A^^) with /3 « 0.5 given by 



Padoan et al. 



(l997| l. 



and a typical value of cr ~ 0.3 = (7in/ln(10) from our results 
for the DIG and warm Hi we obtain M ~ 1.6. |Hill et aLlp008l > 
compared their observed Emission Measure PDFs for the DIG with 
PDFs derived from isothermal MHD turbulence simulations to find 
reasonable agreement for 1.4 < A4 < 2.4, consistent with our 
estimated value. A typical DIG temperature of T ~ 8000 K gives 
Cs ~ 12.5 km and M — 1.6 would then require turbulent 
velocities of u ~ 20kms~^. 
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While the global-disc simulations of |Wada & Normari|j2007^ 
produced lognormal density PDFs, their dispersions are about 4 
times greater than the dispersions we have found. This may be re- 
lated to the average gas densities in their simulations: (riHi) ^ 
lcm~'' compared to e.g. (riHi) ~ 0.3 cm~^ for our total Hi 
sample. 

We may draw the following conclusions from the discussion 
of our results: 

(i) The density PDF of the diffuse ISM cannot be fitted by one 
lognormal. 

(ii) The density PDFs of the diffuse ISM in the disk (|6| < 5°) 
and away from the disk 5°) are lognormal, but the positions 
of their maxima and the dispersions differ 

(iii) Several effects seem to influence the shape of the PDF. An 
increase of the number of clouds/cells along the LOS causes a de- 
crease in the dispersion and a shift of the maximum to higher den- 
sities. On the other, hand, an increase in the average density (or 
decrease in the mean temperature) increases the dispersion as well 
as the density of the maximum. 

The competing effects described in the last conclusion will 
complicate the interpretation of PDFs observed for external galax- 
ies. 

Elmegreen (2002|l and |Wada & Norman| (^007 ) have shown 
that the star formation rate in a galaxy is related to the shape of 
the density PDF if it is lognormal: the dispersion is an important 
parameter in this respect l |Tassis|[2007[ |Elmegreenl|2008^ . A bet- 
ter understanding of the factors that influence the dispersion of the 
lognormal density PDF may be obtained from future simulations. 



5 SUMMARY 

Lognormal density PDFs have been found in recent numerical sim- 
ulations of the ISM - both local, isothermal models (e.g.|Vazquez-| 
ISemadeni & Garcia|200ll|Ostr iker et al. 2001 ; Kowal et al. 2007 1 
and multi-phase, global models |de Avillez & Breitschwerdt)2005( 
|Wada & Norman|2007[ > - and have become an important compo- 
nent of theories of star formation ( [Elmegr een " 2002[ |Tassis|2007[ 
[Elmegreen 2008 1 . To date there has been little observational data 
with which to compare density distributions produced by the simu- 
lations. 

The results reported here provide strong support for the exis- 
tence of a lognormal density PDF in the diffuse (i.e. average den- 
sities of n < 1 cm"'') ionized and neutral components of the ISM. 
In turn, the form of the PDFs is consistent with the small-scale 
structure of the diffuse ISM being controlled by turbulence. Future 
simulations should allow the calibration of the dispersion of the dif- 
fuse gas density PDF in terms of physically interesting parameters, 
such as the number of turbulent cells along the line of sight. 
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